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v  o'  o 
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Force 
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Gravity  constant 
Buoy  drai t 
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Strain  as  a  function  of  (T  -  T  ) 
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Scutum 
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Numerical  factor  =  pC^A^,  (c^  -  cffl)/2 


Numerical  factor  =  pgirdg  /4 

Total  number  of  cable  segments  =  total  number  of  bodies 
Poisson’s  ratio 

Fluid  velocity  relative  to  the  cable  =  c  - 
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Distance  measured  along  the  cable 
Cable  tension 
Drift  velocity 
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Weight  per  unit  length  in  air 


Buoyancy  force  per  unit  length 


Weight  in  water  of  the  bottom  unit 


Direction  normal  to  a  cable  segment 
Horizontal  direction  positive  to  the  right 
Direction  tangential  to  a  cable  segment 
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Angle  with  the  vertical  of  a  chord  connecting  the  top  of  the 
array  with  the  bottom  weight 

p  Fluid  mass  density 

$  Angle  of  the  cable  with  the  horizontal 

Angle  of  the  cable  with  the  vertical 

SUBSCRIPTS 
a  Above 

b  Below 

B  Bottom  unit 

H  Horizontal  direction 

I  Intermediate  body 

i  Irteger  index 

H  Maximum 

m  Minimum 

o  Reference  state 

p  Prescribed  value 
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V  Vertical  direction 

x  x  direction 

y  y  direction 
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ABSTRACT 

An  analysis  is  given  for  the  two-dimensional  steady-state 
behavior  of  extensible  free-floating  cable  systeEC.  The  anal¬ 
ysis  includes  the  differential  equations  of  equilibrium  for  an 
extensible  cable  and  the  iteration  schemes  for  the  equilibrium 
system  drift  velocity  and  surface  buoy  draft.  Based  on  this 
analysis,  a  FORTRAN  IV  program,  FF2E,  has  been  written.  The 
program  allows  the  cable  system -to  have  an  arbitrary  number  of 
different  cable  segments  and  intermediate  bodies.  The  program 
works  well  tei  the  large  majority  of  cases  of  practical  interest. 
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INTRODUCTION 

Until  recently,  steady-state  cable  studies  have  been  concerned 
largely  with  towing  and  single  point  moored  cables.  For  cases  of  a  cable 
towing  a  given  body  or  of  a  single  cable  mooring  a  surface  buoy  with  a 
given  draft,  the  problem  reduces  to  an  initial  value  problem  where  the 
differential  equations  of  equilibrium  of  the  cable  are  integrated  along  the 
cable  for  given  initial  values  at  one  end  of  the  cable.  The  integration 
stops  when  cable  scope  or  the  prescribed  depth  is  reached.  With  the  use 
of  subroutines  for  the  solution  of  differential  equations,  such  as  the 
subroutine  KUTHER  which  uses  the  Kutta-Merson  method,  these  initial  value 
problems  can  be  conveniently  solved  on  modem  day  digital  computers.  At 
NSRDC,  Cuthill1  has  developed  a  program  for  the  two-dimensional  shape  of 
inextensible  towed  and  moored  cables  in  a  uniform  stream.  More  recently, 
Wang2  has  developed  a  program  for  the  three-dimensional  shape  of  an  exten¬ 
sible  cable  system  moored  in  the  presence  of  a  current  profile  which  may 
change  in  magnitude  and  direction  with  depth. 


References  are  listed  on  page  23. 
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Perhaps  the  simplest  case  where  iteration  techniques  are  required 
in  the  steady-state  analysis  of  cable  systems  is  the  case  of  a  cable  of 
given  scope  moored  in  a  given  depth.  In  this  case  the  submergence  of  the 
surface  buoy  must  be  iterated  until  the  cable  reaches  the  given  depth  with 
the  given  cable  scope.  Where  the  current  is  severe  and/or  the  depth  is 
large,  the  buoy  may  be  completely  submerged  beneath  the  ocean  surface. 

The  use  of  iteration  techniques  is  required  in  the  steady-state 
analysis  of  some  modem  applications  of  cable  systems.  Skop  and  O'Hara3 
have  studied  the  problem  of  subsurface  buoys  moored  by  an  array  of  cables. 
The  need  for  iteration  procedures  for  this  problem  can  be  readily  seen  by 
considering  the  case  of  a  single  subsurface  buoy  moored  by  more  than  one 
cable.  It  is  clear  that  the  subsurface  buoy  can  be  maintained  in  equilib¬ 
rium  by  an  infinite  combination  of  the  forces  in  the  cables  just  below  the 
buoy.  Conceptually,  the  forces  at  the  tops  of  these  cables  must  be  con¬ 
tinually  iterated  until  the  lower  ends  of  these  cables  reach  the  required 
anchor  positions.  Instead  of  iterating  the  forces  at  the  top  of  the 
cables,  Skop  and  O'Hara  choose  the  equivalent  method  of  iterating  the 
forces  at  the  lower  end,  i.e.,  at  the  anchor  positions. 

In  the  case  of  a  cable  towed  in  a  circular  path  the  radius  of  the 
motion  of  the  towed  body  must  be  continually  iterated  until  the  boundary 
conditions  at  the  towing  end  are  satisfied.  This  problem  has  been  con¬ 
sidered  by  Huang,4  Choo,5  Skop,0  and  Skop  and  Choo.7 

In  the  case  of  free-floating  cable  systems,  the  problem  considered 
in  the  present  study,  iteration  schemes  are  also  required.  Basically,  the 
drift  velocity  of  the  entire  cable  system  and  the  draft  of  the  surface  buoy 
must  be  continually  varied  until  the  entire  system  is  in  an  equilibrium  of 
forces  in  both  the  horizontal  and  vertical  directions. 

The  present  report  gives  a  description  of  computer  program  FF2E, 
which  is  a  FORTRAN  IV  program  for  the  two-dimensional  steady-state  behavior 
of  extensible  free-floating  cable  systems.  The  iteration  schemes  for  the 
system  drift  velocity  and  buoy  draft  are  described  in  detail.  It  is  shown 
that  the  schemes  work  well  for  the  large  majority  of  cases  of  practical 
interest.  For  some  cases  of  long  cables  in  unusually  severe  current  pro¬ 
files  the  schemes  may  fail  to  converge  due  to  the  extreme  sensitivity  of 
cable  configuration  to  changes  in  drift  velocity  and  buoy  draft.  The 
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differential  equations  of  equilibrium  used  in  the  present  study  are  similar 
to,  but  not  identical  with,  the  equations  used  by  Pode.8  The  differences 
are  simply  due  to  the  difference  in  coordinate  systems  used.  The  program 
allows  the  c^ule  system  to  have  an  arbitrary  number  of  different  cable 
segments  and  intermediate  bodies.  The  surface  buoy  is  taken  to  be  a  circu¬ 
lar  cylinder  in  the  present  study.  However,  by  properly  adjusting  the 
values  for  the  length,  diameter,  and  drag  coefficient,  other  buoy  shapes 
may  be  approximated.  For  cases  where  the  free-floating  system  contains  a 
lower  array  of  hydrophones,  the  program  calculates  the  inclination  and 
curvature  of  the  array. 

In  the  present  study,  the  inclination  of  the  buoy  is  taken  to  be 
vertical.  In  actual  cases,  the  buoy  inclination  is  a  function  of  such 
factors  as  the  center  of  gravity  of  the  buoy,  the  center  of  buoyancy,  the 
point  of  attachment  of  the  cable  to  the  buoy,  and  the  magnitude  and  loca¬ 
tion  of  the  resultant  drag  force  (which  is  dependent  on  the  inclination  of 
the  buoy)  acting  on  the  buoy.  The  inclination  of  the  buoy  must  be  deter¬ 
mined  by  solving  a  transcendental  equation  for  each  net;  assumed  set  of 
values  for  drift  velocity  and  buoy  draft.  This  would  result  in  increased 
complexity  of  program  logic  and  computer  time  requirements .  It  is  felt 
that  the  present  simple  approach  for  the  surface  buoy  will  yield,  in  most 
cases,  reasonable  answers  for  such  system  variables  as  system  drift 
velocity,  cable  configuration,  and  operating  depth  of  the  bottom  unit. 

The  equilibrium  inclination  of  the  surface  buoy  can  then  be  calculated 
from  the  equilibrium  values  of  the  drift  velocity  and  buoy  draft  furnished 
by  the  present  computer  program.  For  those  cases  where  the  effect  of  sur¬ 
face  buoy  inclination  substantially  affects  cable  system  variables,  then 
this  effect  should,  of  course,  be  incorporated  into  the  program. 

The  program  may  also  be  used  to  approximately  solve  cable  payout 
problems.  In  these  applications,  an  object  may  be  kept  stationary  on  the 
ocean  bottom  by  having  an  intermediate  cable  pack  pay  out  cable  as  it  drifts 
along.  The  amount  of  cable  which  is  payed  out  depends  on  the  drift  velocity 
of  the  cable  pack.  The  actual  drift  of  the  cable  pack  is  rather  complex. 

The  weight  of  the  cable  pack  decreases  as  cable  is  being  payed  out,  and  the 
tension  exerted  on  the  cable  pack  by  the  cable  being  payed  out  is  a  func¬ 
tion  of  time.  An  approximate  analysis  of  the  motion  of  the  cable  pack  may 


be  made  by  assigning  average  values  to  the  cable  pack  weight  and  the  verti¬ 
cal  and  horizontal  components  of  tension  exerted  on  the  cable  pack  by  the 
cable  being  payed  out.  The  resulting  cable  system  consisting  of  surface 
float,  connecting  cable,  and  cable  pack  may  then  be  considered  as  a  free- 
floating  system. 

GENERAL  METHOD  OF  SOLUTION 

The  tree-floating  cable  system  considered  in  the  present  study  is 
shown  in  Figure  1.  The  system  consists  of  a  surface  buoy,  an  arbitrary 
number  of  different  cable  segments  and  intermediate  bodies,  and  a  terminal 
weight  at  the  bottom.  As  shown  in  Figure  1,  a  current  profile  which  varies 
with  depth  may  be  prescribed  over  the  entire  cable  system. 

In  the  presence  of  the  given  current  profile,  the  entire  cable  sys¬ 
tem  will  drift  at  a  velocity  Up  and  the  surface  buoy  will  taxe  a  draft  H 
such  that  the  entire  cable  system  is  in  an  equilibrium  of  forces.  The 
principal  problem  of  the  present  study  is  to  determine  these  equilibrium 
values  of  Up  and  H.  Using  the  iteration  schemes  described  later  in  the 
study,  values  for  Up  and  H  are  chosen.  From  tiiese  values  of  UQ  and  H,  the 
current  relative  to  the  cable  system  as  well  as  the  tension  and  angle  of 
the  cable  iust  below  the  buoy  are  determined.  The  differential  equations 
of  equilibrium  may  then  be  integrated  down  the  cable.  At  intermediate 
points  along  the  cable  the  integration  of  the  equations  must  be  interrupted 
to  take  into  account  intermediate  attached  bodies  which  may  be  drogues  or 
hydrophones . 

The  integration  of  the  equations  proceeds  in  the  above  manner  until 
the  bottom  unit  is  reached.  At  this  point,  the  equilibrium  of  the  bottom 
unit  is  checked.  Since  the  conditions  at  the  top  of  the  cable  are  deter¬ 
mined  by  equilibrium  considerations  for  the  surface  buoy,  and  the  cable 
configuration  is  obtained  by  considering  the  equilibrium  of  the  cable,  then 
the  entire  system  is  in  equilibrium  if  the  bottom  lit  is  in  equilibrium. 
Figure  1  shows  that  the  imbalance  of  forces  in  the  horizontal  and  vertical 
directions  may  be  written  as 

eH  =  TH  "  <DB  +  TBH) 
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SJSBS^S 


ev  =  Tv  "  ^WB  +  TBV^ 


(2) 


where  e^  and  e^  are  the  imbalances  in  the  horizontal  and  vertical  direc¬ 
tions,  respectively, 

7^  and  Ty  are  the  horizontal  and  vertical  components  of  the  tension 
just  above  the  bottom  unit,  respectively, 

Dg  and  Wg  are  the  drag  and  weight  in  water  of  the  bottom  unit, 
respectively,  and 

TgH  and  TgV,  in  the  case  of  a  payout  problem,  refer  respectively  to 
the  horizontal  and  vertical  components  of  tension  exerted  by  the 
cable  being  paid  out  on  the  bottom  unit;  for  free-floating  systems, 

TBH  =  TBV  =  °* 

As  mentioned  previously,  the  quantities  Tg^  and  Tg^  are  taken  to  be  con¬ 
stants  in  the  present  study. 

A  number  of  criteria  on  the  magnitude  of  the  errors  may  be  used  to 
determine  when  the  problem  is  considered  to  be  solved.  For  example,  one 
possible  criterion  would  be  to  say  that  if  the  resultant  imbalance  in 
forces  E,  defined  by 


*  V^H 


2  2 

+  ev 


(3) 


is  less  than  a  presciibed  amount  Ep,  the  problem  is  considered  to  be 
solved.3  If  E  >  Ep,  then  new  values  must  be  found  for  Up  and  H.  In  the 
present  study  the  problem  is  considered  to  be  solved  if  the  following  two 
inequalities  are  simultaneously  satisfied 


m 


*eH<  <  el  ' °B  +  TBH‘ 


(4a) 


I  eV‘  <  £2  +  TB\T 


(4b) 


where  and  are  small  positive  constants. 

EQUATIONS  OF  EQUILIBRIUM 

Figure  2  shows  the  coordinate  systems  used  in  the  present  study. 

The  (x,y)  spatial  coordinate  system  is  taken  to  have  fixed  directions,  with 


....  aa  -  »  ----  -  CSiL HtSi* *  -  J  to: 


.  texijLiSu* tfateA  *6**i*. 


the  origin  at  the  undisturbed  ocean  surface.  The  x  axis  is  positive  to 
the  right  and  the  y  axis  is  positive  downwards.  The  cable  differential 
equations  are  written  for  a  coordinate  system  attached  to  the  cable,  called 
the  cable  coordinate  system  (X,Y) .  The  Y  axis  is  directed  along  the  cable. 

CABLE  EQUATIONS 

Differential  Equations 

The  two-dimensional  differential  equations  of  equilibrium  for  a 
flexible  cable  have  been  derived  by  a  number  of  authors;  see,  for  example, 
the  derivation  given  by  Pode.8  For  the  coordinate  systems  used  in  the 
present  study,  the  equations  of  equilibrium  for  the  cable  in  an  arbitrary 
stretched  condition  are 

—  +  G  +  W  sin  <j>  =  0  (5) 

T  +  I  +  W  cos  <p  =  0  (6) 

where  T  is  the  tension, 

s  is  the  stretched  cable  scope, 

W  is  the  weight  per  unit  length  in  water  of  the  stretched  cable,  and 
G,  I  are  the  fluid  drag  forces  per  unit  length  acting  on  the  cable 
in  the  Y  and  X  directions,  respectively. 

The  above  equations  are  similar  to,  though  not  identical  with,  the  equa¬ 
tions  derived  by  Pode.8  The  differences  are  due  to  the  difference  in 
coordinate  systems  used.  The  present  coordinate  system  defines  the  verti¬ 
cal  coordinate  as  being  positive  downwards  whereas  in  Reference  8  it  is 
positive  upwards.  Also,  for  0  <  <J>  £  90  deg,  cable  scope  in  the  present 
study  is  positive  measured  from  the  surface  buoy,  whereas  in  Reference  8 
cable  scope  is  positive  measured  from  the  bottom  weight.  The  present 
coordinate  system  is  more  convenient  for  free-floating  and  moored  cable 
systems  while  the  coordinate  system  used  in  Reference  8  is  more  convenient 
for  towing  cables. 

As  pointed  out  in  Reference  2,  it  is  convenient  to  relate  the  quan¬ 
tities  ds,  G,  I,  and  W  to  quantities  at  some  reference  state  where  all  the 
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cable  parameters  are  known.  This  reference  state  is  taken  to  occur  at 


T  =  T  ,  where  Tq  need  not  be  equal  to  zero.  Summarizing  the  results  of 


Reference  2,  the  relationships  between  ds,  G,  I,  and  W  in  terms  of  these 


quantities  at  the  reference  state  T  =  Tq  are  given  by 


ds  =  ds  (1  +  e) 
o 


I  =  f  (e)  I 
cv  J  o 


G  =  f  (e)  G 
CK  o 


W  =  —22 —  w  [f  (e)]2 
1  +  e  bo  1  cv 


where  subscript  o  denotes  quantities  for  the  reference  state. 


e  is  the  strain  =  K(T  -  Tq) , 


(10a) 


K  is  an  experimentally  measured  tension-strain  function. 


fc^  "  J- 

o 


(10b) 


d  is  the  cable  diameter. 


WaQ  is  the  weight  per  unit  length  in  air  of  the  cable, 


is  the  buoyancy  force  per  unit  length 


WHn  "  P  8 


(10c) 


p  is  the  fluid  density,  and 

2 

g  is  the  gravity  constant  =  32.2  ft/sec  . 

In  the  present  study  the  cable  is  assumed  to  be  linearly  elastic  so  that 


K  and  f  are  given  by 


K  =  Kj  •  (T  -  Tq) 


(11a) 
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fc(e)  -  1  -  pe 


(lib) 


where  K^  is  a  constant  and  p  is  Poisson1 s  ratio.  In  Equation  (9),  W  is 
given  in  terms  of  WaQ  and  W^.  Usually  the  cable  weight  is  specified  in 
terms  of  Wq,  the  weight  per  unit  length  in  water  of  the  cable  at  the  refer¬ 
ence  state.  In  these  cases,  WaQ  is  simply  given  by 


W  =  W  +  W, 
ao  o  bo 


where  is  given  by  Equation  (10c). 

Upon  using  Equations  (7) -(9),  Equations  (5)  and  (6)  become 


-s—  +  G  i  (e)  (1  +  e)  +  sm 
ds  o  cv  J  K  J 
o 


in  {,ao  -  W| 


bo  c 


f  (e)(l+e)  =0  (13) 


a 

%  '(■ 


T  SD  *  Zo  fcte)tl  *  e)  *  cos  *["ao  -  "bo  fc2(e)(1  *  'j]  -  0  <14J 

In  addition  to  the  dependent  variables  T  and  <f>,  other  dependent 
variables  are  of  interest.  These  include  x,  the  horizontal  displacement 
of  the  cable  measured  from  the  point  of  attachment  of  the  cable  to  the 
buoy,  y,  the  vertical  displacement  of  the  cable  measured  from  the  sea 
surface,  and  the  stretched  cable  scope  s.  From  Figure  2,  it  is  seen 
that  the  differential  equations  for  x  and  y  are  given  by 


■jg-  =  -  cos  <j)  (1  +  e) 
o 


■gj-  =  sin  <t>  (1  +  e) 
o 


The  differential  equation  for  s  is  obtained  by  simply  rewriting  Equation 

(7) 


o 


'XVI-X  ^17^' 7 i'*.'-**"*  ^-^,w~Jt  to 

«\  «'  -j>  ' 


Equations  (13) -(17)  along  with  Equation  (10a)  constitute  six  equations 
for  the  six  dependent  variables  T,  <p,  x,  y,  s,  and  e.  Upon  substituting 
Equation  (10a)  into  Equations  (13) -(17),  the  result  is  five  differential 
equations  for  the  five  dependent  variables  T,  <j>,  x,  y,  and  s. 

Fluid  Drag  Forces 

For  bare  cable,  the  normal  drag  force  Iq  is  usually  written  as8 


Io  =  1/2  p  CDdQ  q  sin  <p  |q  sin  <{>| 


'  (18) 


where  CQ  is  the  drag  coefficient,  q  =  c(y)  -  UQ,  and  c  is  the  ocean  current 
which  is  a  function  of  depth. 

There  is  considerable  variation  among  the  forms  proposed  for  the 
tangential  drag  acting  on  bare  cable;  see  the  survey  article  by  Casarella 
and  Parsons.9  In  the  present  study,  the  following  form  proposed  by  Pode8 
is  used 


u  =  ik  — : - 

O  0  jcos  (j)| 

where  f  is  an  empirical  constant  =0.02  for  normal  bare  cable  and 

R  =  1/2  p  Cnd  q I q I .  It  should  be  noted  that  for  the  low  velocities 
o  u  o 

relative  to  free-floating  and  moored  cable  systems,  the  tangential  drag  is 
usually  much  smaller  than  the  other  fluid  and  gravity  forces  occurring  in 
Equations  (13)  and  (14).  Thus,  it  is  not  necessary  to  use  a  complex  form 
for  the  tangential  drag.  Other  forms  for  the  tangential  drag  can,  of 
course,  be  conveniently  incorporated  into  the  formulation. 

EQUATIONS  FOR  INTERMEDIATE  BODIES 

The  integration  of  the  preceding  equations  must  be  stopped  when  an 
intermediate  body  such  as  a  drogue  or  hydrophone  is  reached.  The  tension 
and  angle  just  below  the  body,  T^  and  <f>^ ,  must  be  such  that  the  body  is 
in  an  equilibrium  of  forces.  These  forces  are  the  tensions  in  the  cable 
above  and  below  the  body  and  the  fluid  and  gravity  forces  acting  on  the 
body.  Upon  solving  the  two  equations  of  equilibrium  for  the  body  the 
following  expressions  for  T^  and  are  obtained 
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T  sin  <p 
a  Ta 


T  cos  6 
a  Ya 


(20b) 


where  T  and  <J>  are  the  tension  and  angle  of  the  cable  just  above  the  body. 

Si  Si 

Fjx  is  the  drag  force  acting  on  the  body,  and 
is  the  weight  in  water  of  the  body. 

The  drag  force  FJx  is  given  by 


FIX  =  i/'2  P  cdAI  q|qI 


(21) 


where  CpAj  is  the  drag  area  of  the  body  for  horizontal  flow. 


ITERATION  SCHEMES 

For  chosen  values  of  UQ  and  H,  the  integration  c4'  the  equations  in 
the  preceding  section  proceeds  until  the  bottom  unit  reached.  At  this 
point,  the  equilibrium  of  the  bottom  unit  is  checked.  As  pointed  out 
previously,  if  the  bottom  unit  is  in  equilibrium,  then  the  entire  cable 
system  is  in  equilibrium.  Practically  speaking,  if  the  errors  e^  and  ey 
are  small  enough  to  satisfy  the  given  error  criterion,  the  problem  is 
considered  solved.  If  the  criterion  is  not  satisfied,  then  new  values 
must  be  found  for  Up  and/or  H. 

Computer  program  FF2E  contains  a  preliminary  scheme  and  two  prin¬ 
cipal  schemes  for  iterating  values  of  Up  and  H.  They  are  as  follows, 
in  the  order  that  they  appear  in  the  program: 

(a)  Preliminary  scheme  for  Up. 

(b)  Simultaneous  binary  scheme. 

(c)  Staggered  binary  scheme. 

These  schemes  are  described  below. 
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PRELIMINARY  SCHEME  FOR  Up 

In  this  scheme,  only  the  value  of  Up  is  varied.  The  value  of  H  is 
fixed  and  is  determined  by  the  buoyancy  required  to  support  the  weight  of 
the  entire  system  in  water,  in  the  absence  of  any  current. 

The  scheme  searches  for  a  value  of  Up  which  gives  a  value  of  e^ 
less  than  some  small  prescribed  quantity  e..  The  scheme  does  this  by 
restricting  Up  to  lie  between  successively  closer  barriers.  For  eH  greater 
(less)  than  zero,  AUp  is  chosen  so  that  the  new  value  of  Up  lies  halfway 
between  the  present  value  of  Up  and  the  most  recent  preceding  value  of  Up 
which  is  greater  (less)  than  the  present  value  of  Up.  Basically,  the 
changes  in  Up  are  based  on  the  sign  given  in  Equation  (28),  which  is 
derived  in  the  following  section. 

The  principal  purpose  of  the  present  scheme  is  to  provide  an  improved 
initial  value  of  Up  for  the  simultaneous  binary  scheme.  It  has  been  found 
that  with  this  initial  value  of  Up,  the  convergence  of  the  simultaneous 
binary  scheme  is  improved.  Also,  for  certain  cases  of  low  ocean  currents, 
the  sready-state  draft  differs  very  little  from  the  value  of  H  fixed  in  the 
present  scheme,  i.e.,  the  value  of  e„  associated  with  the  final  value  of 
Up  is  also  small.  Thus,  in  these  cases,  the  present  scheme  presents  a 
direct  solution  to  the  problem.  Finally,  as  will  be  seen  later,  the  present 
scheme  essentially  is  the  first  cycle  of  the  staggered  binary  scheme. 

SIMULTANEOUS  BINARY  SCHEME 

This  scheme  changes  the  values  of  H  and  Up  simultaneously. 

Consider  first  the  forces  in  the  horizontal  direction.  The  term 
(Th  -  Dg)  appearing  in  the  formula  for  eH,  given  in  Equation  (1),  repre¬ 
sents  the  horizontal  component  of  the  drag  forces  acting  on  the  entire 
cable  system.  It  should  be  noted  that  (TH  -  Dg)  is  a  complex  quantity 
which  depends  on  the  current  relative  to  the  cable  and  the  inclination  of 
the  cable.  An  approximate  measure  for  (TH  -  Dg)  is  given  by 

th  -  db  ”  1/2  0  Vt  [1/2  (cm  -  V2  -  1/2  ‘UD  -  c/]  (22) 
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II 


where  C^  =  IC^  ♦  £  CQi  lQ±  dQ.  * 


i?!  CDAIi)  » 


C^Ag  is  the  drag  area  of  the  surface  buoy  and  any  underwater  packages 
mounted  on  the  buoy, 

N  is  the  total  number  of  cable  segments  =  total  number  of  bodies, 

S)i*  ^oi*  an<*  ^oi  are  resPectively  t^ie  drag  coefficient,  reference 
length,  and  refexence  diameter  of  the  ith  cable  segment, 

CpAj^  is  the  drag  area  of  the  ith  body, 

c^  is  the  maximum  value  of  the  current,  and 


c  is  the  minimum  value  of  the  current, 
m 


II  can  be  seen  from  the  above  definition  for  C^Ap  that  it  is  a  measure 
of  the  drag  area  of  the  entire  cable  system.  Noting  that  is  simply 
a  constant  in  Equation  (1),  it  can  be  seen  that  e^  depends  directly  on 
CTH  -  Dg)*  In  particular,  it  can  be  seen  that  if  eH  is  greater  (less)  than 
zero,  the  change  in  (T^  -  Dg) ,  A(TH  -  D_) ,  should  be  less  (greater)  than 
zero.  Thus,  the  error  efJ  may  be  related  to  A(TH  -  Dg)  as  follows 

eH  *  '  «TH  -  V  '  '  1/2  0  CD*T  [S  '  V'-  %>  -  %  -  c„>  %] 

-  “d  P  Vt  1/2  CcM  -  c»>  -  mJ\  <23> 


where  M,,  =  p  1/2  (cM  -  cm 
then  given  by 


(cu  -  c  ) .  An  expression  for  AU_  in  terms  of  e„  is 
Mm  u  H 


Consider  now  the  vertical  forces.  By  noting  that  the  quantities 
Wg  and  Tgy  are  constants  in  the  formula  for  ey  given  in  Equation  (2),  it 
can  be  seen  that  e^  depends  directly  on  Ty  In  particular,  it  can  be  seen 
that  if  ev  is  greater  (less)  than  zero,  the  change  in  Ty  ATy  should  be 
less  (greater)  than  zero.  Provided  that  changes  in  the  vertical  component 
of  tension  just  below  the  buoy  are  propagated  without  a  reversal  in  sign 


-.»<*.'W«*  AS:**-*'  y*:-'*1'.*  Sv^3* 


l| 

£ 


E 


«-. ...  -c^%*>. 


L 


B  / 


to  the  bottom  of  the  cable,  e^f  may  be  related  to  changes  in  the  buoyancy  of 
the  surface  buoy  Bg.  The  buoyancy  force  Bg  is  given  by 


ird„ 


BS  =  pg  ~ 


H  = 


*v* 


(25) 


where  dg  is  the  diameter  of  the  surface  buoy  and  My  =  pg  Trdg 2^4.  From  the 
preceding  discussion,  ev  may  be  related  to  ABg,  the  change  in  Bg,  as 


follows 


ev  =  -  ATV  *  -  ABg  =  -  My  AH 


(26) 


An  expression  for  AH  in  terms  of  ey  is  then  given  by 


AH  -v 


«V 


(27) 


A  number  of  proportionality  factors  have  been  tested  for  Equations 
(24)  and  (27) .  No  factors  were  found  which  cause  the  scheme  to  converge 
to  the  equilibrium  values  of  H  and  Up  for  all  cases.  The  optimum  defini¬ 


tions  for  AUp  and  AH  were  found  to  be  of  the  following  form 


AUd  = 


6  •  SH  h)!  eH 


E  • 


”u 


(28) 


AH  = 


*  !ey|  ey 


E  •  M, 


(29) 


where  6  is  a  positive  proportionality  quantity,  0  <  6  £  1  and  is  a 
positive  quantity.  As  discussed  in  the  following  paragraph,  the  quantity 
6  serves  the  purpose  of  reducing  the  changes  in  Up  and  H  if  these  changes 
lead  to  a  resultant  error  which  is  larger  than  the  preceding  resultant 


error.  The  quantity  6^  is  initially  set  equal  to  1.  If  the  change.;  in  e^ 


occur  too  sivrffly  (fast),  then  is  increased  (decreased).  The  ratios 
je^J/E  and  |ey|/E  which  occur  in  Equations  (28)  and  (29),  respectively, 
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5^* 


basically  tell  the  scherte  to  concentrate  on  changing  Up(H)  xf  the  predomi¬ 
nant  error  is  e^Oy) .  It  was  found  that  the  presence  of  these  radios 
increases  the  range  of  convergence  of  the  scheme. 

The  scheme  starts  by  choosing  the  initial  values  for  H  and  Up  to 
be  the  final  values  of  these  quantities  from  the  preliminary  scheme  and 
calculating  the  cable  configuration  from  these  initir.l  values.  The 
quantities  AUp  and  AH  are  computed  from  the  resulting  error  values  e^,  e H, 
and  H  by  using  an  initial  value  of  6  equal  to  1.  These  values  of  AUp  and 
AH  are  added  respectively  to  Up  and  H,  resulting  in  new  values  of  H  and 
Op,  H'  and  Up'.  The  cable  configuration  is  again  calculated  using  the  new 
values  H'  and  Up*.  New  values  of  error  are  obtained:  e^' ,  ev',  and  E’. 

If  E*  <  E,  then  6  remains  at  its  initial  value,  and  the  new  values  H'  and 
Up'  are  retained.  The  quantities  AUQ  and  AH  are  obtained  from  e^* ,  ev' , 
and  E'.  if,  however,  E*  >  E,  then  the  new  values  e^’ ,  e^’,  E‘,  H',  and 
Up*  are  rejected.  The  quantities  AUp  and  AH  are  obtained  from  the  old 
values  of  eH,  Cy,  and  E  using  a  smaller  value  of  5.  The  quantity  6  is 
continually  reduced  until: 

(a)  E'  <  E, 

(b)  e^'  simultaneously  satisfies  the  following  two  conditions: 

(eH,/eH5  >  1  3,1(1 

{eH'i  >  lev'^  or 

(c)  6  is  lower  than  a  prescribed  lower  bound.  This  lower  bound 
is  set  equal  to  0.2  for  the  first  ten  iterations  of  the  scheme  and  O.OS 
for  subsequent  iterations. 

In  any  of  the  above  three  cases,  5  is  reset  equal  to  its  initial  va^ue  of 
1,  and  AUp  and  AH  are  computed  by  using  the  new  error  values  e^’ ,  ev', 
and  E'.  The  resulting  values  of  AUp  and  AH  are  added  to  Up'  and  H' , 
respectively.  The  scheme  continues  in  the  above  manner  until  the  errors 
eH  and  ey  are  sufficiently  small  so  that  the  given  error  criterion  is 
satisfied. 

In  one  of  the  earlier  schemes  tested,  the  quantity  &  was  allowed  to 
decrease  indefinitely  until  E’  <  E.  In  other  words,  options  (b)  and  fcj 
of  the  preceding  paragraph  were  not  present.  Furthermore,  6  was  never 
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reset  equal  to  its  original  value,  but  either  remained  the  same  or  decreased 
(according  to  whether  E'  was  less  than  or  greater  than  E)  as  the  iteration 
proceeded.  This  scheme  was  similar  to  one  proposed  by  Skop  and  O'Hara, 3 
who  dealt  with  the  problem  of  subsurface  floats  moored  by  an  array  of  cables. 
It  was  soon  found  that  this  scheme  did  not  converge  in  many  cases.  It  was 
found  that  at  certain  steps  of  the  iteration,  the  new  error  E*  could  never 
be  made  smaller  than  E,  no  matter  how  small  5  was  made.  It  was  found  that 
this  was  due  to  interaction  effects  whereby  changes  in  H  have  large  effects 
on  e^  and/or  changes  in  Up  have  large  effects  on  ey.  It  should  be  noted 
that  interaction  effects  are  not  represented  in  Equations  (28)  and  (29) . 

It  was  found  that  these  interaction  effects  arise  in  cases  where 
the  current  differential  between  the  top  and  bottom  ot  the  cable  system  is 
fairly  large.  In  these  cases,  the  normal  fluid  drag  forces  acting  on  the 
cable,  given  in  Equation  (18) ,  became  significantly  large.  The  interaction 
effects  arise  largely  due  to  the  strong  dependence  of  the  drag  forces  on 
cable  geometry  as  shown  in  Equation  (18).  Consider,  for  example,  the  case 
where  only  Up  or  H  is  changed.  For  a  given  change,  the  geometry  of  the 
entire  cable  may  be  changed  significantly,  and  the  resulting  changes  in  the 
normal  drag  forces  may  have  significant  components  in  both  directions. 

Skop  and  O'Hara3  prove  that  their  scheme  converges  for  moored  cable 
arrays.  It  should  be  noted,  however,  that  the  proof  depends  on  the  assump¬ 
tion  that  all  of  the  forces  acting  on  the  cable  are  unaffected  by  changes 
in  the  geometry  of  the  cable.  Thus,  their  work  should  be  viewed  with  some 
caution,  particularly  for  cases  of  high  ocean  currents  where  the  geometry 
dependent  fluid  drag  forces  become  significant. 

The  results  of  Reference  3  do  suggest,  however,  that  the  convergence 
of  the  present  scheme  may  be  improved  by  restricting  the  cable  shape  to 
lie  between  certain  bounds.  An  auxiliary  scheme  has  been  incorporated  into 
the  program  which  changes  Up  by  small  amounts  whenever  the  angle  $  exceeds 
or  becomes  less  than  certain  prescribed  bounds.  This  auxiliary  scheme  has 
been  found  to  improve  the  convergence  of  the  main  scheme.  It  should  be 
noted  that  the  bounds  on  <{>  can  neither  be  too  close  nor  too  far  apart.  If 
the  bounds  on  $  are  too  far  apart,  then  it  is  clear  that  the  scheme  loses 
its  effectiveness.  On  the  other  hand,  excessively  close  bounds  may  mean 
many  changes  of  Up,  thus  increasing  computer  time  requirements.  Also,  the 


equilibrium  cable  configuration  may  lie  outside  the  bounds  if  the/  are  too 
restrictive.  It  has  been  found  that  the  optima  lower  and  upper  bounds  for 
$  are  approximately  0  degrees  and  125  degrees,  respectively. 

As  aentioned  previously,  this  scheae  does  not  converge  for  all  cases 
The  limits  of  convergence  see®  to  depend  aost  strongly  on  the  current  dif¬ 
ferential  between  the  top  and  bottom  of  the  cable  systea,  and  on  the  length 
of  the  cable.  Based  on  many  hundreds  of  runs  already  made  with  the  program 
it  seems  that  the  scheme  converges  for  the  large  majority  of  cases  where 
the  current  differential  is  less  than  3.0  knots.  In  these  cases,  less  than 
approximately  twenty  iterations  are  usually  required  to  obtain  values  of 
H  and  Up  which  yield  values  of  the  resultant  error  E  less  than,  say,  0.05 
pounds.  For  some  cases  of  high  current  differentials,  the  scheae  cannot 
minimize  E  beyond  a  certain  amount,  ranging  from  a  few  tenths  of  a  pound 
to  several  pounds.  In  these  cases,  the  computer  program  switches  to  the 
staggered  binary  scheme. 

STAGGERED  BINARY  SCHEME 

Description  of  Scheae 

The  preliminary  scheme  constitutes  the  first  cycle  of  the  staggered 
binary  scheme.  Basically,  the  scheae  changes  H  and  Up  one  at  a  time,  thus 
avoiding  the  interaction  effects  aentioned  above.  The  changes  in  H  and 
Up  are  made  in  accordance  with  the  signs  indicated  in  Equations  (29)  and 
(28),  respectively.  First,  an  initial  value  of  H  is  chosen  to  be  the 
value  required  to  support  the  weight  of  the  entire  system  in  water,  in  the 
absence  of  any  current.  Then,  keeping  the  value  of  H  constant,  the  scheme 
searches  for  a  value  of  Up  which  gives  a  value  of  e^  less  than  e^.  The 
procedure  for  doing  this  has  already  been  discussed  in  the  section  dealing 
with  the  preliminary  scheae. 

When  a  value  of  Up  has  been  found  which  gives  rise  to  <  ej»  the 
iteration  for  Up  stops.  The  scheae  calculates  a  new  value  for  H  based  on 
the  value  of  e^,  the  value  of  ey  associated  with  e^  <  e^,  as  follows 
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ah  = 


"5V  evo 


«v 


where  6y  is  a  positive  number  and  My  is  defined  in  Equation  (25). 

It  may  be  noted  that  Equation  (30)  differ?  from  Equation  (29)  in 
that  the  quantities  |ey|/£  and  6  which  occur  in  Equation  (29)  have  been 
replaced  by  1  and  6y,  respectively.  For  the  case  of  eyg,  where  <  e^, 
the  ratio  |ey|/E  is  nearly  equal  to  1.  The  quantity  6y  is  initially  chosen 
to  be  0.6.  Should  the  convergence  to  the  equilibrium  value  of  H  be 
excessively  slow,  5y  is  increased.  Similarly,  5y  is  decreased  if  the 


changes  in  H  cause  excessively  large  changes  in  eyp.  It  should  be  noted 


that  H  is  also  restricted  to  be  between  successively  closer  barriers.  For 
eyg  greater  (less)  than  zero,  the  new  value  of  H  is  restricted  to  be 
greater  (less)  than  the  last  preceding  value  of  H  which  gives  rise  to  an 


error  eVQ  less  (greater)  than  zero.  Should  the  value  of  AH  calculated  in 


Equation  (30)  lead  to  a  new  value  of  H  which  goes  beyond  the  preceding 
barrier,  the  new  value  of  H  is  adjusted  so  that  it  lies  between  the  pre¬ 


sent  value  of  H  which  gives  rise  to  the  givsn  value  of  eyo  and  the  barrier. 


Keeping  this  new  value  of.  It  constant,  the  scheme  again  searches  for  a  value 
of  Up  which  gives  rise  tc  eH  <  e^.  The  srneme  proceeds  in  the  above  manner 
until  the  error  quantities  e^  and  ey  are  -.mall  enough  to  satisfy  the  given 
error  criterion 


Discussion  of  Convergence 


This  scheme  converges  over  a  wider  range  of  cable  and  current  param¬ 
eters  than  the  simultaneous  binary  scheme.  It  is,  however,  slower,  usually 
requiring  at  least  sixty  iterations  to  arrive  at  suitably  small  values  of 
eu  and  e„.  Based  on  many  hundreds  of  computer  runs  made  thus  far.  it  seems 
that  the  scheme  converges  for  nearly  all  cases  where  the  surface  current 
is  less  than  approximately  5.0  knots.  It  should  be  noted  that  surface 
currents  seldom  exceed  5  knots  (see,  for  example.  Reference  10). 

For  some  cases  of  long  cables  floating  in  the  presence  of  surface 
currents  in  excess  of  5  knots,  the  scheme  may  fail  to  converge.  In  these 
cases,  where  the  geometry  dependent,  fluid  drag  forces  dominate,  the  scheme 
cannot  minimize  the  error  E  beyond  a  certain  amount  ranging  from  a  few 
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tenths  of  a  pound  to  several  pounds.  The  computer  results  for  these  cases 
show  that  the  cable  configuration  and  the  error  quantities  e^  and  are 

very  sensitive  to  changes  in  H  and  Up.  As  an  example,  consider  the  partic¬ 
ular  case  of  a  2,000-ft  cable  floating  in  the  presence  of  a  surface  current 
of  8  knots.  For  this  case,  the  computer  results  showed  that  at  H  =  0.86318 
ft,  a  change  in  H  in  the  sixth  or  higher  decimal  place  produced  a  change  in 
e^Q  from  -6.2  lb  to  l.l  lb.  This  large  change  of  e^  from  negative  to 
positive  values  is  typical  of  the  behavior  of  the  scheme  when  it  fails  to 
converge. 

The  preceding  results  were  obtained  for  the  case  when  the  differ¬ 
ential  Equations  (13) -(17)  were  integrated  to  an  accuracy  of  0.01  percent. 
In  view  of  the  extreme  sensitivity  of  e^0  to  changes  in  H,  it  was  decided 
to  investigate  this  particular  case  further  by  increasing  the  accuracy  of 
the  integration  of  the  differential  equations.  With  accuracies  up  to  and 
including  0.0001  percent,  the  large  changes  of  with  changes  in  H  in 
the  sixth  or  higher  decimal  place  around  e^Q  =  0  continued  to  occur.  The 
values  of  H  giving  rise  to  these  large  changes  in  eVQ  changed  somewhat  with 
the  accuracies.  For  example,  with  accuracies  of  0.001  percent  and  0.0001 
percent,  the  large  changes  in  eVQ  from  negative  to  positive  values  occurred 
respectively  at  H  =  0.86199  ft  and  H  =  0.86066  ft.  The  investigation  had 
to  be  terminated  when  an  attempt  was  made  to  integrate  the  differential 
equations  with  an  accuracy  of  0.00001  percent.  Computer  time  requirements 
became  prohibitively  large. 

The  preceding  discussion  shows  the  nature  of  the  physical  and/or 
numerical  instability  for  cases  of  unusually  severe  current  profiles  and 
sufficiently  long  cables  where  the  cable  geometry  dependent  fluid  drag 
forces  dominate.  In  order  to  thoroughly  investigate  the  nature  of  this 
instability,  and  to  determine  the  solution,  computer  time  requirements 
beyond  the  scope  of  this  study  would  be  required.  The  preceding  results 
do  tentatively  indicate  that  the  equilibrium  value  of  H  is  localized  to 
a  small  neighborhood  by  the  scheme.  In  most  cases  where  the  scheme  fails 
to  converge,  the  minimum  resultant  error  does  not  exceed  a  few  pounds.  It 
is  expected  that  in  most  engineering  applications  involving  cibles  a  few 
thousand  feet  or  longer,  an  error  of  one  or  two  pounds  can  be  tolerated. 


DISCUSSION  OF  UNIQUENESS 


The  uniqueness  of  the  numerical  solutions  has  not  been  investigated 
in  the  present  study.  One  way  to  study  this  question  is  to  incorporate  a 
grid  scheme  into  the  computer  program  whereby  H  and  U^  are  each  assigned  a 
given  number  of  values  over  prescribed  intervals.  The  configuration  of  the 
cable  and  the  resultant  error  E  are  then  calculated  for  each  possible  pair 
of  values  of  H  and  Up.  If  E  is  equal  to  zero  for  two  or  more  different 
pairs  of  values,  then  the  solutions  are  nonunique.  Based  on  the  sensitiv¬ 
ity  of  E  to  changes  in  H  and  UQ,  observed  from  the  many  runs  using  program 
FF2E,  it  seems  that  in  most  cases  a  minimum  of  50  to  100  values  must  be 
prescribed  for  both  H  and  U^,  around  their  respective  equilibrium  values 
predicted  by  program  FF2E,  in  order  to  discover  the  existence  of  any  addi¬ 
tional  solutions  which  may  exist.  This  means  that  for  a  given  case,  the 
configuration  of  the  cable  must  be  calculated  2,500  to  10,000  times.  The 
execution  time  on  the  CDC  6700  for  a  single  calculation  of  the  configuration 
of  the  cable  is  of  the  order  of  one  second.  Thus,  several  hours  of  computer 
time  probably  would  be  required  in  order  to  investigate  uniqueness  for  one 
case.  A  systematic  investigation  would  require  the  consideration  of  a  large 
number  of  cases.  The  computer  times  required  were  beyond  the  scope  of  the 
present  study. 

FORTRAN  IV  COMPUTER  PROGRAM 

Based  on  the  preceding  analysis,  a  FORTRAN  IV  computer  progrf-m, 

FF2E,  has  been  written  for  the  system  shown  in  Figure  1.  The  program  is 
compatible  with  the  CDC  6700  computer  currently  being  used  at  NSRDC.  The 
program  consists  of  a  main  program  and  four  subroutines.  The  main  program 
reads  in  the  parameters  for  the  surface  buoy,  cable  segments,  intermediate 
bodies,  and  bottom  weight,  as  well  as  the  wind  speed  and  current  profile. 

The  main  program  also  prints  out  the  input  data. 

It  is  worthwhile  to  make  several  comments  about  the  data  which  are 
read  in.  The  program  reads  in  data  for  a  cylindrical  surface  buoy.  How¬ 
ever,  by  reading  in  suitable  values  for  the  length,  diameter,  and  drag 
coefficient,  other  surface  buoy  shapes  may  be  approximated.  Also,  as 
previously  mentioned,  the  program  assumes  that  the  cable  is  linearly 


elastic  by  reading  in,  for  each,  cable  segment,  the  constant  quantity  K^, 
given  in  Equation  (11a).  It  is  felt  that  this  approximation  is  sufficiently 
accurate  for  the  large  majority  of  engineering  applications.  If  one 
wishes  to  accurately  model  nonlinear  tension-strain  behavior,  the  subroutine 
ELAS  which  is  Dresent  in  program  MR3E2  may  be  conveniently  incorporated  into 
the  present  program.  It  should  be  noted  that  in  the  above  program  the 
differential  equations  are  integrated  only  once.  As  pointed  out  earlier  in 
the  section  on  iteration  schemes,  the  differential  equations  may  have  to  be 
integrated  100  times  or  more  in  the  present  program  before  a  solution  is 
achieved.  Thus,  the  incorporation  of  the  subroutine  ELAS  will  increase 
computer  time  requirements  for  the  present  program.  Finally,  it  should  be 
noted  that  in  some  free-floating  cable  systems,  the  lower  intermediate 
bodies  form  an  array  of  hydrophones.  In  this  case,  the  reader  simply  reads 
in  the  number  of  lower  intermediate  bodies,  exclusive  of  the  bottom  weight, 
which  form  the  array.  As  pointed  out  later  in  the  present  section,  the 
computer  program  prints  out  additional  information  about  the  configuration 
of  the  array. 

After  reading  in  the  input  data,  the  program  searches  for  the. 
equilibrium  values  of  Up  and  H.  This  is  done  by  the  iteration  schemes, 
described  previously,  contained  in  the  subroutine  STEADY.  For  given 
values  of  H  and  Up,  the  differential  Equations  (13)-(17)  are  integrated 
down  the  cable.  The  differential  equations  are  defined  in  the  subroutine 
DAUX  and  integrated  by  the  subroutine  KUTMER,  which  uses  the  Kutta-Merson 
method  for  solving  systems  of  first-order  differential  equations.  The 
differential  equations  are  integrated  to  an  accuracy  of  0.01  percent.  By 
simply  changing  one  or  two  statements  in  the  program,  any  accuracy  may  be 
prescribed.  As  mentioned  earlier,  computer  time  requirements  increase  with 
increasing  accuracy.  The  proper  value  of  the  current  is  given  by  the  sub¬ 
routine  CUR  which  linearly  interpolates  the  input  current  profile  points. 
The  calculation  of  Equations  (20a) -(21)  in  the  case  of  intermediate  bodies 
is  performed  by  the  subroutine  STEADY. 

When  the  bottom  of  the  cable  is  reached,  the  errors  e^  and  e^  are 
calculated  and  are  tested  to  see  if  they  satisfy  the  inequalities  given 
in  Equations  (4a)  and  (4b).  If  both  inequalities  are  not  simultaneously 
satisfied,  new  values  must  be  found  for  H  and/or  Up.  Before  doing  sc. 
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the  subroutine  STEADY  prints  out  the  iteration  number,  the  values  of  Up, 

H,  ev,  e^,  Dg,  6  (used  in  Equations  (28)  and  (29)),  and  the  tensions  in 
the  cable  just  below  the  surface  bucy  and  just  above  the  bottom  weight. 

If  both  of  the  inequalities  (4a)  and  (4b)  are  simultaneously  satis¬ 
fied,  the  program  stops  iterating  for  further  values  of  H  and  UQ  and  goes 
on  to  print  the  following  information: 

1.  System  drift  velocity  Up  and  buoy  draft  H. 

2.  Vertical  and  horizontal  components  of  the  tension  in  the  cable 
just  below  the  surface  buoy. 

3.  Variables  at  points  along  the  cable,  as  functions  of  the  refer¬ 
ence  cable  scope  sQ  measured  from  the  surface  buoy: 

s,  stretched  cable  scope, 

x,  horizontal  displacement  relative  to  buoy  bottom, 

y,  vertical  displacement  relative  to  ocean  surface, 

T,  tension  in  the  cable, 

<p,  angle  measured  from  the  horizontal,  and 

<j>v,  angle  measured  from  the  vertical. 

4.  Fluid  and  gravity  forces  Dg  and  Wg  acting  on  the  bottom  unit. 

5.  Vertical  and  horizontal  components  of  the  tension  just  above 
the  bottom  unit. 

In  the  event  that  there  are  two  or  more  lower  hydrophones,  exclusive 
of  the  bottom  weight,  the  program  prints  out  0,  the  angle  with  the  vertical 
of  a  chord  connecting  the  top  of  the  array  with  the  bottom  weight,  and  k, 
the  maximum  perpendi cu] ar  distance  between  the  chord  and  array.  The  quan¬ 
tities  0  and  k  are  shown  in  Figure  3.  It  should  be  noted  that  the  coor¬ 
dinate  system  (xr,  yr)  shown  in  Figure  3  is  defined  such  that  the  bottom 
weight  is  the  origin  of  the  coordinate  system. 

After  printing  out  the  above  information,  the  program  goes  on  to  a 
new  case. 
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